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Abstract 

We study an extended Gross-Neveu model with Nf quark flavors and with 
an additional SU(2) (global color) degree of freedom of the quarks. The four 
fermion interaction in the color channel is mediated by a random color matrix 
with fixed strength. The effective potential for the quark condensate and the 
scalar correlation function is investigated in the large iVj-limit. Quark anti- 
quark thresholds are absent in the scalar correlation function implying that 
the decay of the scalar meson in free quarks is avoided. The transition from 
the low energy to the high energy perturbative regime is found to be smooth. 
Perturbation theory provides a good description of Green's functions at high 
energies, although there is no quark liberation. 



* Supported in part by DFG under contracts La-932/1 — 1/2. 
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1 Introduction 



One of the most challenging problems nowadays in hadron physics is to understand 
why mesons and baryons do not decay into free quarks (quark confinement), al- 
though these are the fundamental constituents of Quantumchromodynamics (QCD), 
the right theory of strong interactions. Lattice simulations verify the quark confine- 
ment from a theoretical point of view |], j2]. More and more evidence arose in the 
recent past, that the quark confinement is due to a condensation of monopoles which 
subsequently expels electric flux from the vacuum ||. This provides a linear rising 
potential between static quarks implying confinement. A recent success by Seilberg 
and Witten |4], H shows that monopole condensation is indeed the mechanism of 
quark confinement in certain super-symmetric Yang-Mills theories. 
Despite the increasing knowledge of the low energy properties of QCD, hadron 
physics is still only feasible by effective quark models f|, [7J. They successfully 
describe the mechanism of the spontaneous breakdown of chiral symmetry and the 
physics of the light mesons, but inherently suffer from a non- confinement of the 
quarks J7|. Recently, an effective quark model was proposed which incorporates 
the quark confinement in a sense that quark anti-quark thresholds are absent in 
(mesonic) Green's functions ||. The main ingredient in this model is an effec- 
tive low energy four quark interaction mediated by a random gluonic background 
field. The ground state properties of this model were investigated under extreme 
conditions, i.e. high temperature and/or density. The results were found to be in 
agreement with the phenomenological expectations. 

In order to study further the impact of random background fields on the quark 
confinement, we here study an extended version of the Gross- Neveu model ||. As in 
QCD, the effective (running) coupling strength of this model becomes small at high 
energies (asymptotic freedom). Here, we will illustrate a mechanism which serves as 
a possible answer to the question how perturbation theory at high energies provides 
a good description of Green's functions, although there is no quark liberation. 

2 Ground state properties 

2.1 Model description and effective potential 

The model under investigations is described by the Euclidean generating functional 
in two dimensions for mesonic Green's functions, i.e. 

m = {l,jv q v- q exp{-/ A [L - HWM] }) Q , (1) 
L = q(x)(0 + im)q(x) + ^-{q{x)q{x)) 2 + ^—qr a q(x) qr h q(x) G a \ (2) 
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where m is the current quark mass. In the limit m = 0, the theory is invariant 
under global chiral rotations. The quarks fields q(x) transform under a global O(Nf) 
flavor and a global SU(2) (color) symmetry, g and G ab are dimensionless coupling 
constants in the color singlet and the color triplet channel respectively. An average 
of all orientations O of the background field G ab , transforming as G' = T GO with 
O being a 3 x 3 orthogonal matrix, is understood in (TfJ) to restore global SU (2) color 
symmetry. From || we expect that these randomness leads to the confinement of 
quarks. 

In order to explore ground state properties of the model, it is convenient to study the 
effective potential U of the (Euclidean) condensate, since a non-vanishing value sig- 
nals the spontaneous breakdown of chiral symmetry. The effective action T([gqq] c ) 
of the fermion condensate is obtained by a Legendre transformation of Z[j] with 
respect to j, i.e. 

f 5 Z [j] 

r([0??]c) := -Z[j] + J d 2 x [gqq]c(x)j(x) , [g qq\ c {x) := jjr~r . (3) 

The effective potential U results from F by restricting [g qq] c {x) to constant values, 
i.e. U([gqq] c ) = r([gqq] c )/V, where V is the Euclidean four volume. The poten- 
tial U is determined up to an unphysical offset, which is usually chosen in order 
the potential U to be zero at vanishing condensate [gqq] c - The value of the quark 
condensate at which the effective potential attaches its global minimum is the vac- 
uum expectation value (g qq) . The corresponding minimal value of U is the vacuum 
energy density. Serveral minima correspond to different phases of the model. The 
state with the lowest energy density (value of U at the local minimum) forms the 
true vacuum. 

2.2 Solution of the model in the large iV^-limit 

Although an exact solution of the model might be feasible with the powerful tech- 
niques developed in the context of the Gross- Neveu model |10| , it is sufficient for our 
purposes to study the model (|l|) in leading order of a l/iV/-expansion. In the context 
of two-dimensional models, the 1/iV-expansion was cast into doubt at least for mod- 
els with a continuous symmetry like the chiral Gross-Neveu model, since it predicts 
a spontaneous breakdown of the symmetry. The occurrence of Goldstone bosons 
in two dimensions causes infra-red problems implying that the symmetry cannot 
be broken spontaneously []rT| . However, it was subsequently shown by Witten that 
the symmetry is "almost" spontaneously broken, and that the large iV-expansion 
provides a rather good guide to the ground state properties |T2| . 
The 1 / iVy-expansion is most easily derived by rewriting the theory (|I]) by means of 
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auxiliary fields , i.e. 

Z\j\ = (hx J Vq Vq X?£ X?£ a exp |- J d 2 x [L M - j(x) gq(x)q(x)} }^ , (4) 

L M = g(x)(^ + ^ + zSV a )g(x) + ^(S -m) 2 + ^S^G- 1 )^ 6 . (5) 

Integrating out the quarks in (^) one obtains an equivalent theory formulated in the 
fields S , S a only. The non-perturbative approach applied here is an expansion with 
respect to 1/Nf for a fixed background field G ab and an average over all orientations 
O of G ab once the Green's function has been calculated. For this purpose, we 
decompose the auxiliary fields into a mean-field part and fluctuations, i.e. 

E (x) = M + a (x) , £ a (x) = iM a + a a (x) , (6) 

where the orientation of the mean-field M a depends on the actual choice of G ab . 
The imaginary unit in front of M a in was introduced for convenience. Later 
it will turn out that the ground state favors an imaginary constituent quark mass 
(corresponding to real M a ) in the color triplet channel. This will be the main obser- 
vation implying the absence of quark anti-quark thresholds in mesonic correlation 
functions. 

We will evaluate the functional integral in (^) in a semi-classical expansion around 
the mean field values. (It was recently pointed out that mean-field theory might 
generically provide a good description of the ground state of fermionic theories even 
at finite temperatures |13|)- As in the case of the standard Gross-Neveu model 0, 



fluctuations <To, o a are suppressed by a factor 1/Nf implying that we might confine 
ourselves to the mean-field level (as long as ground state properties are addressed). 
For constant external sources j, it is sufficient to consider constant mean field values. 
A straightforward calculation yields 

1 -Z[j] = -^-Trln{i$ + i{M + iM a T a )) (7) 



VN f V 

+ i-(M - m - ig 3 ) 2 - l -M a {G- l ) ab M b , 

where the trace extends over internal indices as well as over space time. The master 
fields Mo and M a are solutions of the gap equations 

dm _ dz[j] _ 

It is easy to check that the right hand three equations reduce to an eigenvector 
equation to determine M a , i.e. 



G ab M b = g c M a , M : v 7 .U".\/" . (9) 
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and an equation to calculate M. There are three eigenvalues of the matrix G ab . 
With g c we define the eigenvalue of G ab which subsequently provides the phase with 
the lowest action. Using (pj), the trace term in (0) can be reduced to 



[ A du In (V + 2u{Ml - M 2 ) + (M 2 + M 2 ) 2 
4n Jo v 



7r — 2arctan 



A. 



A 2 - A 2 



const. 



A- In 41 
A 4 



(10) 



where A± = Mq±M 2 and A is the 0(4)-invariant cutoff to regularize the momentum 
integration. 

In order to renormalize the model, one observes that the logarithmic divergence in 
(jin) is proportional to A_ and can be therefore absorbed by a redefinition of the 
coupling constants, i.e. 



1 , A 2 1 1 
- In — + - = — , 

TT IF 9 9R 



1 , A 2 1 1 
-In— + - = — , 

TT /i 2 g c g R 



(11) 



In order to completely renormalize the theory, we also define the renormalized mass 
by 

m tub 

— = — . 12 

g gR 

Field renormalization is not requested at the present stage of approximation. The 
/3-function is in either case of the coupling constants (3(gii) = —2/iig R and signals 
asymptotic freedom. The renormalized coupling constants g^, g c R decrease by in- 
creasing momentum scale /i. Note that in coincidence the bare coupling constants 
vanish in the infinite cutoff limit. This implies that the terms m 2 /g, gj, mj approach 
zero for A — » oo. 

We are now ready to perform the Legendre transformation (|3|). The calculation is 
straightforward. The final result for the effective potential of the quark condensate 
in the large A^-limit is 



1-U 



2g R 

2vr 



gR 



2g c R 



-M 



Ml - M 2 / M 2 + M 2 



2ty 



In 



(13) 



7T 



M 2 

2 arctan ■ 



M 2 



2M M 



where M is directly related to the quark condensate by 



M 



(14) 
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and M satisfies the gap-equation 

7T %e „ f M 2 + M 2 M /tt M 2 -M 2 

■I + I In —2— + — H - - arctan H = . (15) 



g%(ji) [ /i 2 M \2 2M M 



The equations 013,14,15) are one of our main results. We have obtained a renormal- 



ization group invariant result, since a change of the renormalization point [i can be 

absorbed by a redefinition of the renormalized coupling constants. 

In order to remove the superficial dependence of the effective potential U on the 

subtraction point //, we introduce two renormalization group invariant scales Sq, s c 

by 

11s 2 lis 2 

+ -ln^| = 0, — l -- + -ln% = 0. (16) 



g R (ji) tt /i 2 ' g%(fj) tt /i 2 

With the help of (|i~6"l), the effective potential becomes 

1 tt m R^ M l M 2 + M 2 




f M 2 + M 2 M /tt M 2 - .,, 

M \ In — ^ + — - - arctan —2 U = . 1 

1 s 2 M I 2 2M M 1 r 




where 



Thereby equation ([TaD was used to simplify the expression for the effective potential. 
The two free parameters gn, g c R of the original theory are removed in favor of the 
two renormalization group invariant scales sq, s c (dimensional transmutation). One 
of these scales, e.g. so, can be used as fundamental unit to express dimensionful 
quantities. The remaining free parameter is the ratio of these scales, or lnsg/s 2 . 
Equation (|T8|) is numerically solved for given value of M . The potential U is 
subtracted in order to yield zero at M Q = 0. The numerical result for the effective 
potential U in (|17]) is shown in figure 1. Note that M = is always a solution of 
([TE|). In this case, the effective potential U ( |TTD coincides with that of the standard 
Gross-Neveu model (solid line in figure 1). Note, however, that for sufficiently small 
values of M Q a solution with M/0 exists. This solution has lower effective potential 
(dashed line in figure 1). The global minimum of this branch therefore constitutes 
the vacuum. The effective potential U is smooth at the transition from the phase 
with M / to the phase M — 0, although M varies rapidly (dot-dashed line in 
figure 1). 

Ground state properties are described by the values of the quark condensate at 
which the effective potential attaches its global minimum. A non-zero value of M 
corresponds to a non- vanishing quark condensate (see (13)), whereas M controls the 
absence of quark anti-quark thresholds in Green's functions as we will see below. In 
order the effective potential to have a minimum, M and M must satisfy 

m R f, M 2 + M 2 M fn M 2 -M 2 \] , , 

7r — — = M I hi - arctan — 9 } , 19 

g R \ si M ^2 2M M I J 1 ; 
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in addition to flllf) . Figure 2 shows the lowest action solution. If the overall scale 
is chosen to be Sq, the only parameter is lns^/s^. One observes that if s c exceeds a 
certain value, i.e. 

s c > 2.71828. . . s , (20) 

a phase with a non-vanishing value of M occurs. For values of s c less than this 
critical value, M is zero and the phase coincides with that of the standard Gross- 
Neveu model. In section |3], we will discuss the remarkable properties of this phase. 



2.3 Quark propagator and analytic structure 

Here, we will discuss the Euclidean quark propagator in the phase characterized by 
M ^ 0. In momentum space, it is 

S{k) = ft + i{M + iM"r«) ' (21) 
Introducing the eigenvectors |±) of the matrix M a T a , i.e. 

M a r a \±) = ±M\±) , (22) 



the quark propagator fl2T|) decomposes into two parts with conjugate complex masses, 
i.e. 

1 1 

s{k) = l+) p + Z (m + iM) (+l + H jrm^M) H ' (23) 

The hermiticity and unitarity properties of a theory with a pair of particles with 



conjugate complex masses were extensively studied in several models |14|, [TfJ, [T? . 

One finds evidence that such a theory does not necessarily violate unitarity. 

In order to evade the triviality problem [Il5j of the scalar 4 -theory, a Higgs model 



with two scalar degrees of freedom with complex conjugate masses was proposed |]TB 
A consistent field theory is obtained by quantizing the scalar modes in a Hilbert 
space with indefinite metric. The model can be described by an Euclidean path 
integral, although the Minkowskian functional integral might not exist due to the 
presence of the so-called ghost particles fTI| . It is argued that microscopic acausality 



effects remain undetectable using realistic wave packets, and that the S-matrix of 



the model is unitary fT6 



Recently, the connection of a confining potential and the occurrence of complex sin- 
gularities in the propagator of the fundamental degrees of freedom was investigated 
in QED3 |7|. 

From the literature, we conclude that a model containing particles with conjugate 
complex masses does not necessarily violate the phenomenological requirements for 
causality and unitarity of the S-matrix. In contrast, such a model might be a 
candidate for a low-energy theory of QCD describing the confinement of quarks. 
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In this paper, we will not further question the analytic properties of the extended 
version of the Gross-Neveu model studied here. Although this model lacks an im- 
mediate application to hadron physics, we hope it is simple enough to allow for 
a clarification of the phenomenology of particles with conjugate complex masses. 
These investigations are left to future work. 



3 The scalar correlation function 

In order to focus on the physical impact of the imaginary constituent quark mass in 
the color triplet channel, we here concentrate on the scalar correlation function, i.e. 

A(p 2 ) = J d'x e-*- (gq q ( x ) gqq(Q)) = J d 4 x li=* ■ (24) 

Note that Mq and therefore gqq is the renormalization group invariant and the 
physical quantity. Unphysical quark anti-quark thresholds manifest themselves as 
imaginary part of the correlation function. In the following, we will establish the 
absence of an imaginary part of A(p 2 ) in the phase with M^O. 
In order to obtain A(p 2 ) in leading order of the l/iV/-expansion, we expand the 
Lagrangian @ up to second order in the fluctuations in @, i.e. 

Z\j] « (in J Va Va a exp{-S( 2 )} ^ , (25) 
S (2) = J-0~ 2 {lvm°(p>(-p) + \cr a mf(p>"(-p) (26) 



ia a (p)K a (p 2 )a(-p) - i^tj(p)j(- p ) -iN f j(p)a(-p) 



where j(p) is the Fourier transform of the external source j(x) and 

N f r d 2 k 
~~g~ ~ J (2tt^ 



AT, f Ail- 

n°(p 2 ) = - tr {S(k + p)S(k)} , (27) 



Uf(p 2 ) = NfiG- 1 )^ - J -0- 2 tr {r a S(k + P yS(k)} , (28) 

K a ( P 2 ) = ~i J-^- 2 tr{T a S(k + p)S(k)} . (29) 

The quark propagator S(k), in momentum space, is given by (21). The calculation 
of the functions ( pT| - p9]) is straightforward and closely follows [§]. It turns out that 
the quantities (p7|) and (|29|) can be expressed in terms of two functions H (p 2 ), 
H v (p 2 ), i.e. 

AT AT A/fa 

n°(p 2 ) = ^ H (p 2 ) , K«{p 2 ) = H v (p 2 ) . (30) 

7T 7T M 



s 



One finds 



ffo P 2 = - + - da In °- + 2 , 31 

g 2 Jo A 4 

H v (p 2 ) = — da arccos . , 
J o ,Jw 2 + AM 2 Mq 

where 

w = a(l- a)p 2 + Ml - M 2 . (32) 

Whereas H v is finite, H Q needs renormalization. Equation (pi]) , which was employed 
to renormalize the effective potential, absorbs the cutoff dependence in fl3~T|) as ex- 
pected. Removing the superficial dependence of Ho(p 2 ) on the subtraction point [i 
with help of the renormalization group invariant scale sq in flTB"!), one finally obtains 
the finite and explicitly renormalization group invariant result 



H (P 2 ) = i; L dot ln^— - A °- + 2 , (33) 



1 ri w 2 + AM 2 M 2 
- / da In 4 

2 Jo Sq 

We also find that M a is an eigenvector of the polarization matrix 11"", i.e. 



Equipped with the results ( PD| ) and (P4|), it is now easy to integrate out the meson 



fields do and a a in ( p5[) and to calculate (|24j). The result of this procedure only 
depends on M 2 (besides p 2 and Ml). This implies that average over all orientations 
of the background field G ab in (|25|) is trivial, since M is a singlet. The final result 
for the scalar correlation function A(p 2 ) in leading order of the l/iVy-expansion is 
therefore 



_KN f 
H o(P 2 ) + H {pi)+\nsl/sl 



We are now going to study the occurrence of an imaginary part of the scalar cor- 
relation function signaling a quark anti-quark threshold. For this purpose, we first 
perform the analytic continuation of the Euclidean correlation function to Minkowski 



space. This continuation follows the lines of JT4J. The correlation function in 
Mikowski-space HS.m{p 2 m) as a function of the Minkowski four-momentum pu is 
related to the Euclidean counterpart by 

iA M (p 2 M ) = A(-p 2 M ) . (36) 

It is therefore sufficient to study the Euclidean correlation functions at negative 
momentum squared. 
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In order to illustrate the disappearance of the quark anti-quark threshold for M^O, 
we first study its occurrence for M = 0. In this case, our model describes the scalar 
correlation function in the usual Gross-Neveu model with a constituent quark 
mass Mq. The term of interest is the integrand of Ho(p 2 ) in which essentially 
becomes Iuw/sq. This implies that whenever w becomes negative, the function Hq, 
(|33|), acquires an imaginary part, which subsequently describes the decay of the 
scalar meson in a free quark anti-quark pair. In order for w to become negative, the 
Euclidean momentum p 2 must satisfy, 

-P 2 < 4M 2 , (37) 

implying that the quark-anti-quark threshold occurs at a (Minkowskian) momentum 
Pm = 2Mo, which is the familiar result. For M = the function H v (p 2 ) does no 
harm, since it is identically zero. 

In contrast, for M^O the functions H / V (p 2 ) are real for all momentum p 2 . No quark 
anti-quark threshold occurs at all. This is our main observation. Figure 3 shows 
the scalar correlation function as function of p 2 . At negative values of p 2 , a peak 
occurs near the would-be threshold position. This peak is enhanced, if the system 
is driven towards the deconfinement phase by a choice of the parameter lns^/s^. 
The correlation function smoothly approaches the perturbative result (dashed line 
in figure 3) at large momentum transfer. Our model therefore represents an example 
the Green's functions of which can be accurately calculated in perturbation theory 
at high momentum, although there is no quark liberation. 



4 Conclusions 

An extended version of the Gross-Neveu model [[J for Nf quark flavors was studied. 
In addition to the 0(iV/) -flavor symmetry, the model possesses a global SU(2) (say 
color) symmetry. The interaction in the color triplet channel possesses a random 
orientation in color space. The effective potential for the quark condensate was cal- 
culated in the large Nf limit. Two phases were found: in one phase, the constituent 
quark mass in the color triplet channel is zero (M = 0), and a non- vanishing dynam- 
ical generated mass M in the color singlet channel occurs. This phase corresponds 
to the vacuum phase of the standard Gross-Neveu model. However, a second phase 
constitutes the vacuum, since this phase has lower vacuum energy density. This non- 
trivial phase is characterized by an imaginary constituent quark mass in the color 
triplet channel (M ^ 0). The impact of an imaginary mass on causality of Green's 



functions and unitarity of the S-matrix is discussed in the literature [14], [L6|, \T7j . It 
is argued that these models do not contradict the phenomenological requirements. 
In order to clarify on this issue in the context of the model presented here, further 
investigations are needed, which are left to future work. 
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Rather than to investigate these theoretical aspects, we here focus on the remark- 
able physical consequences of the phase with M/fl. For this purpose, the scalar 
correlation function was studied. In the Gross-Neveu phase (M = 0) of the model, 
an imaginary part of the correlator corresponding to the quark anti-quark thresh- 
old was found. In contrast, the correlation function of the phase with M ^ is 
real, and the quark anti-quark threshold is absent. The occurrence of an imaginary 
constituent quark mass in the color triplet channel might therefore serve as a sig- 
nal for quark confinement (which is defined as the absence of quark thresholds in 
Green's functions). At large momentum transfer, the scalar correlation function can 
be accurately calculated by perturbation theory. The model provides example that 
asymptotic freedom is perfectly compatible with the non-liberation of quarks. 
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Figure captions 



Figure 1: The effective potential U and the color triplet constituent quark mass M 
as function of the quark condensate M in the chiral limit itir = for two values of 
lnsg/s^; all quantities are in units of s . 

Figure 2: The color-singlet and color-triplet constituent quark masses M and M in 
the chiral limit jjir = 0. 

Figure 3: The scalar correlation function as function of the Euclidean momentum 
transfer p 2 . 
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